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Abstract 



We study the properties of the right-handed sneutrino and its viability as a 
WIMP dark matter candidate in an extended version of the NMSSM in which a 
right-handed neutrino superfield is included with a coupling to the singlet Higgs 
in order to provide non-vanishing Majorana neutrino masses. We perform a 
systematic study of the parameter space, including LEP constraints and exper- 
imental bounds on low-energy observables. We investigate the conditions under 
which the right-handed sneutrino has the correct relic abundance and the dom- 
inant annihilation channels. Next we calculate the theoretical predictions for 
the sneutrino-proton elastic scattering cross section and compare it with present 
and future experimental sensitivities. We find that sneutrinos with a mass in 
the range of 5-200 GeV can reproduce the observed dark matter relic density 
without being excluded by direct dark matter searches and for natural values 
of the input parameters. Interestingly, the predicted scattering cross section is 
generally within the reach of future experiments. Finally, we comment on the 
possible implications for collider physics. 



1 Introduction 



A number of astrophysical and cosmological observations have provided substantial 
evidence supporting the existence of dark matter, an abundant component of our Uni- 
verse made of some new, yet undiscovered, particles. A generic weakly interacting 
massive particle (WIMP) is one of the most promising and attractive candidates for 
cold dark matter in our Universe [1], since it would be thermally produced in sufficient 
amount to account for the observed dark matter density. Interestingly, WIMPs appear 
in many well motivated extensions of the standard model providing new physics at the 
TeV scale, such as supersymmetry (SUSY). 

In supersymmetric models a discrete symmetry (R-parity) is often imposed in order 
to forbid lepton and baryon violating processes which could lead, for instance, to rapid 
proton decay. A remarkable phenomenological consequence is that supersymmetric 
particles can only be produced or annihilate in pairs, thus rendering the lightest SUSY 
particle (LSP) stable. If it is neutral and weakly interacting, the LSP can therefore 
be an excellent WIMP candidate for dark matter. Two particles exist with these 
properties in the minimal supersymmetric standard model (MSSM), namely the lightest 
neutralino [2, 3] and the sneutrino [4]. The former is the mixture of neutral gauginos 
and Higgsinos, and one of the most popular and well studied dark matter candidates [5]. 
In contrast, the (left-handed) sneutrino in the MSSM turns out not to be viable as 
dark matter. Being a left-handed field, the sneutrino has a sizable coupling with the 
Z boson, which entails a too large annihilation cross section and therefore a too small 
relic abundance. Furthermore, its scattering cross section off nuclei, also mediated by 
Z boson exchange, is so large that direct detection experiments for dark matter would 
have already observed it [6]. It should be noted that the inclusion of a lepton number 
violating operator can reduce the detection rate [7]. A re-analysis of the sneutrino LSP 
in the MSSM in light of an updated data can be found in Ref. [8]. 

There are motivations to consider extensions of the MSSM. One of them is the obser- 
vation of neutrino oscillation phenomena, suggested by the solar and atmospheric neu- 
trino anomalies and confirmed by various experiments [9] , which imply non- vanishing 
(albeit very small) neutrino masses. However, in the MSSM as well as the SM, neutrinos 
are massless and the tiny neutrino mass is usually explained by a see-saw mechanism in 
which right-handed Majorana neutrinos, N, are introduced [10] with a large Majorana 
mass (usually taken to be around the grand unification scale, but which can also be of 
order of the electroweak scale or even lower). 
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In supersymmetric models, at same time, one has to introduce right-handed sneu- 
trinos, which would then mix with the left-handed sneutrinos. Interestingly, a mixed 
left- and right-handed sneutrino has a reduced coupling to the Z boson thereby al- 
leviating the abovementioned problems and making it viable as dark matter [11, 12]. 
However, a significant left-right mixture is only possible in some particular supersym- 
metry breaking schemes with a large trilinear term [11, 13]. Such a mechanism is not 
available in the standard supergravity mediated supersymmetry breaking, where trilin- 
ear terms are proportional to the small neutrino Yukawa couplings. Recently, another 
realization of large mixing was pointed out [14] by abandoning the canonical see-saw 
formula [10] for neutrino masses. Another possibility is having a pure right-handed 
sneutrino [15, 16, 17, 18, 19, 20, 21, 22, 23]. These cannot be thermal relics, since 
their coupling to ordinary matter is extremely reduced by the neutrino Yukawa cou- 
pling [15, 16, 17, 18], unless they are somehow coupled to the observable sector, for 
example via an extension of the gauge [19, 20] or Higgs [21, 22, 23] sectors 1 . 

The other motivation to extend the MSSM is the so-called "/z problem" [25]. The 
superpotential in the MSSM contains a bilinear term \1H\H2 which gives supersymmet- 
ric Higgs/Higgsino mass terms. Successful radiative electroweak symmetry breaking 
(REWSB) [26] requires 11 of the order of the electroweak scale. The next-to minimal su- 
persymmetric standard model (NMSSM) addresses this problem in a very natural way. 
A new singlet superfield, S, is introduced [27] and the bilinear term (forbidden by a 
global Z 3 symmetry) is promoted to a trilinear coupling \SHiH 2 . Once REWSB takes 
place the singlet field acquires a vacuum expectation value (VEV), v s = (S), thereby 
generating an effective 11 parameter, 11 = \v s . The NMSSM also alleviates the "little hi- 
erarchy problem" of the Higgs sector in the MSSM [28]. Due to the significant changes 
in the Higgs and neutralino sectors, the NMSSM also has an attractive phenomenology 
[29] and interesting consequences for neutralino dark matter [30, 31, 32]. Notice finally 
that although the Z 3 symmetry of the NMSSM may give rise to a cosmological domain 
wall problem, this can be avoided with the inclusion of non-renormalisable operators 
[33, 34]. 

Motivated by above two issues, in Ref. [22] an extension of the MSSM was proposed 
in which two new singlet superfields were included, as in Refs. [35, 21]. An extra 
singlet superfield S addresses the fx problem in the same way as in the NMSSM and 
provides extra Higgs and neutralino states, while an extra singlet superfield N accounts 
for right-handed neutrino and sneutrino states. This model possesses two interesting 

1 Recently, non-LSP right-handed sneutrino dark matter model was also proposed [24]. 
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new consequences. The right-handed neutrino mass is generated spontaneously (or 
dynamically) with the electroweak symmetry breaking. Due to the non-vanishing VEV 
of the singlet Higgs, an effective Majorana mass for the right-handed neutrino appears 
through the new coupling SNN which is of order of the electroweak scale, in the 
same way as the effective \i term is generated in the NMSSM. Moreover, the singlet 
S, which couples to the MSSM Higgs sector, provides electroweak scale interactions of 
the right-handed sneutrino with the rest of the MSSM fields. Thus, the purely right- 
handed sneutrino LSP has the properties of a WIMP. In Ref. [22] it was shown that 
right-handed sneutrinos in this construction can not only be thermally produced in 
sufficient amount to account for the dark matter in the Universe but also that their 
elastic scattering cross section with nuclei is large enough to detect them through Higgs 
exchange processes. 

In this work we perform a systematic analysis of the properties of the right-handed 
sneutrino as a dark matter candidate in the model of Ref. [22]. In Section 2 the theoret- 
ical details of this construction are introduced, explaining in detail the properties of the 
sneutrino spectrum. In Section 3 the thermal production of right-handed sneutrinos is 
analysed for several examples in the parameter space. The right-handed sneutrino relic 
abundance is calculated and the importance of the different annihilation channels is 
investigated in detail. In Section 4 the direct detection properties of the right-handed 
sneutrino are studied, comparing the theoretical predictions of its elastic scattering 
cross section off nuclei with present and future experimental sensitivities. In Section 5 
we briefly comment on the possible implications of this model for collider physics, de- 
scribing the characteristic signals that could be expected. In Section 6 we summarize 
our conclusions. The technical details are relegated to the Appendices. In Appendix A 
we include the Feynman rules for the new diagrams of the model, and in Appendix B we 
detail the calculation of the sneutrino relic abundance and give the explicit expressions 
for the different annihilation channels. 

2 The Model 

The model consists of an extended version of the NMSSM with a new singlet right- 
handed neutrino superfield N. The superpotential is therefore given by 

W = Wnmssm + X N SNN + y N L ■ H 2 N, (2.1) 
WWssm = Y U H 2 ■ Qu + Y d H\ • Qd + Y e E x ■ Le - XSHi • H 2 + j^S 3 , (2-2) 
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where flavour indices are omitted and the dot denotes the SU(2)l antisymmetric prod- 
uct. As in the NMSSM, a global Z 3 symmetry is imposed for each superfield, so that 
there are no supersymmetric mass terms in the superpotential. Note that the term 
NNN and SSN are gauge invariant but not consistent with R-parity and thus are not 
included. Notice also that N does not have a bare Majorana mass but acquires a mass 
through the non- vanishing singlet Higgs VEV, v s . 

The supersymmetric scalar potential for squarks, sleptons, Higgses and the right- 
handed sneutrino, N, is given as V — Vp + Vd with 

V F = lY^u + Y^df + lY^Qf + lY^Q^ + lY^e + yN^N^ + lY^Ll 2 
+ \Y d Qd + y N Le - \SH 2 \ 2 + \Y u Qu - \SH l + y N LN\ 2 
+ | - XH^ + kS 2 + Aat7V 2 | 2 + \2X N SN + y N LH 2 \ 2 , (2.3) 

and 

2/^1 i i i i i \ ^ 

V D = ^l^l—H^Hl-H. + Q^-Q + ^—u + d^-d + U—L + i^] 

n 2 / a _a a a \ 2 

+|EK ffi+fl} y ff2+(5t i (5+i V • (2 - 4) 

The soft SUSY breaking terms are 

->C sca i armass = m 2 ^\Q\ 2 + ml\u\ 2 + m||dj 2 + m\\L\ 2 + m 2 -\e\ 2 

+m 2 Hi \H 1 \ 2 + m 2 H2 \H 2 \ 2 + m 2 \S\ 2 + m%\N\, (2.5) 

where the new soft scalar masses and m s are included, and 

-^A-terms = (A u Y u H 2 Qu + A d Y d H x Q~d + A e Y e H x L~e + H.c.) 

+ {-\A X SH X H 2 + ^A K ,S 3 + H.c.^ 

+ [\ n A Xn SN 2 + y N A yN LH 2 N + H.c.) , (2.6) 

which contains the new trilinear soft terms A\ N and A VN . The sum of the supersym- 
metric and soft SUSY breaking terms give the total scalar potential. 



2.1 Neutrino mass 

As stated above, in this construction, right-handed neutrino masses are generated by 
the non-vanishing VEV of the singlet Higgs as 

M N = 2X N v s , (2.7) 
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being therefore of order of the electroweak scale. Then, in order to reproduce the small 
masses of the left-handed neutrinos, which are given as 



m 



M N ' 



(2.8) 



the low scale seesaw mechanism implies small Yukawa couplings of (9(10~ 6 ) or less. 
Here, v\^ 2 = (-^1,2) denote the VEV of the Higgs doublet. To reproduce light neutrino 
masses and mixing for neutrino oscillation data we would need to introduce the genera- 
tion structure in the right-handed neutrino sector. However, as we will see, these small 
neutrino Yukawa couplings are completely irrelevant for dark matter physics. Hence, 
for simplicity, we consider one generation case, but one may regard that the considered 
sneutrino corresponds to the lightest one among multi-generations. 

2.2 Sneutrino masses 



The sneutrino mass matrix can be read from the quadratic terms with respect to L 
and N 

V(L,N) C \y N H 2 N\ 2 + \2\ N SN\ 2 + \- XSHi+yNLNl 2 
+ | - A#i# 2 + kS 2 + X N N 2 \ 2 + D - term 

+m\\L\ 2 + m%\N\ + (\ n A Xn SN 2 + y N A m LH 2 N + H.c.) . (2.9) 
Decomposing the left-handed sneutrino and right-handed sneutrino iV as 

(2.10) 

y z. 

the sneutrino quadratic term can be written as 

/ *L1 \ 



*>L = ■fypLX + M>L2), N = -L(jV! + iN 2 ), 



-(z> Li , Ni, z> L2 , N 2 )M 2 Sncutlino 
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(2.12) 



Here, we defined 



m L i = m| + |t/Art> 2 | 2 + D — term, 
m LR = yiv(-A^i) f + y N A N v 2 , 

m LR = yNV 2 (-\v s )\ 

m RR = m 2 ^ + \2\ N v s \ 2 + \y N v 2 \ 2 , 

m RR = {A Xn v s + (kv 2 - A^i^) 1 ") . (2.13) 

If these are real, in other words no CP violation, the real part and imaginary part of 
sneutrinos do not mix and its mass matrix (2.11) is simplified as 

l ll m lR + m l 



Eq. (2.11) = Wtfjl ™- + m - 



+ -{y L2 N 2 ) ( m LR + m LR | | | (2 14) 




Note that the mixing between left-handed and right-handed sneutrinos is induced by 
m 2 LR and ttl 2 l ^, both of which are proportional to the neutrino Yukawa coupling y^. 
As mentioned in the previous subsection, yjy must be as small as 0(lCr( 6-7 )) to gen- 
erate sub-eV left-handed neutrino mass via low scale seesaw mechanism according to 
Eq. (2.8). If we rewrite sneutrinos in the mass eigenstates, 

fii = N* >L u L + N£N, (2.15) 

by rotating with an unitary matrix N u , then, unless the diagonal elements of the mass 
matrix (2.14) are extremely degenerated, we find 

Thus, the mixing between left-handed and right-handed sneutrinos is also of (9(1CT 6 — 
1CT 7 )) and therefore negligible. For all practical purposes in this paper, we can ignore 
all off-diagonal elements in Eq. (2.14) and regard sneutrino mass eigenstates as pure 
left- or right-handed fields. One may see that m RR splits the masses of N 1 and N 2 . In 
particular, N 2 is heavier than Ni for m 2 RR < (and viceversa), and this is the situation 
we will consider throughout the rest of this work 2 . 

2 We only do this for simplicity. The imaginary component of the right-handed sneutrino is as good 
a dark matter candidate as the real component. In fact they only differ in the expression for the 
annihilation into a pair of right-handed neutrinos. 
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Figure 1: Sneutrino spectrum and right-handed neutrino mass, M N , as a function of the 
sneutrino soft mass for A\ N = -500 GeV and A^ = 0.01, 0.05. The NMSSM 
parameters have been chosen as A = 0.1, k = 0.05, tan/3 = 5, and ji = 200 GeV. 



Let us briefly illustrate these properties of the sneutrino spectrum. As an explicit 
example, we have taken [i = 200 GeV and A = 0.1 (which implies v s = 2000 GeV), k = 
0.05, together with tan/3 = 5, A\ = 400 GeV, and A K = 0. We fix A Xn = -500 GeV 
and consider two cases, with A^ = 0.01 and 0.1. From Eq. (2.7) the right-handed 
neutrino mass can be calculated to be Mjy = 40 GeV, and 400 GeV, respectively. The 
sneutrino masses are then an increasing function of m^. These features are displayed in 
Fig. 1, where the resulting sneutrino masses are plotted, together with the right-handed 
neutrino mass, as a function of m^. 

As commented above, the mass-splitting of the two sneutrino mass eigenstates 
is dictated by the m 2 RR term in Eq. (2.13). An increase in \n implies not only an 
enhancement of in Eq. (2.7), but also leads to a larger mass difference between Aq 
and N 2 through the increase in This is clearly evidenced in the two panels of 

Fig. 1. Notice that although the neutrino Majorana mass, Mn, always contributes to 
the sneutrino mass through the Tn RR - term in Eq. (2.13), the negativeness of the m RR 
term allows sneutrinos lighter than right-handed neutrinos. 

Variations in k, and A\ N have no impact on the right-handed neutrino mass, but 
affect the mass splitting of sneutrino states through their contribution to the m 2 RR 
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0.05 0.1 0.05 0.1 



Figure 2: (m^, \ N ) plane for the same example as in Fig. 1 but with A x = 250 GeV (left) 
and 500 GeV (right). The trajectories with a fixed sneutrino mass are indicated by means 
of dashed lines. The dotted line represents the trajectory along which the sneutrino mass 
equals the right-handed neutrino mass (above that line, > M N ). 



term in Eq. (2.13). Depending on the relative signs of the various terms in this ex- 
pressions, variations in k and A\ N can enhance or decrease the lightest sneutrino mass, 
respectively. 

In order to illustrate the effect of the A\ N parameter, we have represented in Fig. 2 
the trajectories in the (m^, Ajv) plane with a fixed sneutrino mass for A\ N = —250 GeV 
and —500 GeV. For a given value of the sneutrino mass a larger \A\ N \ allows a wider 
range of values of A at. We will later use this flexibility in order to look for regions in 
the parameter space with the correct relic abundance. 

Variations in A and k also alter the allowed range of values for A at given a fixed 
sneutrino mass, since they affect the m RR term. In general, an increase (decrease) in |A| 
(\k\) leads to a larger |m|j R | and therefore to a larger range in \ N , with an effect that 
mimics the increase in \A\ N \ discussed above 3 . Notice also that, for a given \i = Xv s , 

3 Remembcr in this sense that we are concentrating on the case m 2 RR < 0, which leads to the real 
component of the sneutrino to be the lightest particle. In the opposite case, m RR > 0, when the 
lightest component is the imaginary one, the range in A at would be enlarged for a decrease (increase) 
in |A| (\k\), again due to the enhancement of \m RR \. 
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an increase in | A | entails a reduction of the right-handed neutrino mass through the 
decrease in v s . This makes it easier to obtain sneutrinos with a mass larger than the 
right-handed neutrino, which might be welcome in order to obtain the correct relic 
abundance as we will see in the next section. 

For completeness, we display in Fig. 3 the contours with a constant sneutrino mass 
in the (A, k) plane for a specific choice of Ajv, m^, and A\ N . As |A| (\k\) increases 
(decreases), the splitting between the two right-handed sneutrino states (real and 
imaginary components) decreases and at some point they become degenerate. This 
is indicated by means of a dashed line in the plot. For our choice of parameters the 
real component of the right-handed sneutrino is the light state below that line and 
the imaginary component becomes lighter in the region above. Notice that (depending 
on the values of the input parameters) the sneutrino can become tachyonic (e.g., for 
small values of k or small values of A). Obviously, the tachyonic regions are larger for 
a smaller sneutrino soft mass. 

Finally, despite the importance of variations of tan /3 in the NMSSM phenomenol- 
ogy, they only affect the term V\V 2 in Eq. (2.13), which is proportional to tan/3/(l + 
tan 2 j3) and have therefore little impact on the resulting right-handed sneutrino spec- 
trum. 



2.3 Vacuum and Higgs potential 

So far, we have assumed that (N) = holds in the vacuum. However, the right- 
handed sneutrino may have a non vanishing VEV, depending on the specific choice of 
parameters [35] . We will check here the validity of that assumption and the implications 
for the relevant parameters of the model. 

The part of the scalar potential related with Higgses and sneutrinos is given by 

V = \y N H 2 N\ 2 + \XSH 2 \ 2 + \XSH 1 \ 2 +\- XH^ + kS 2 + X N N 2 \ 2 + \2X N SN\ 2 
+V D (H l7 H 2 ) 

+m 2 Hi \H,\ 2 + m 2 H2 \H 2 \ 2 + m 2 s \S\ 2 + m%\N\ 2 

+ (-XAxSHtfz + ^kA k S 3 + X n A Xn SN 2 + H.c.) . (2.17) 

As showed in the previous subsection, the left-handed and right-handed sneutrinos 
are almost decoupled from each other. Thus, for this discussion, only right-handed 
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Figure 3: Trajectories for constant sneutrino mass in the (A, k) plane the same case de- 
scribed in Fig. 1 and the choices A^ = 0.07, = 100, 150 GeV and A\ N = 250, 
500 GeV. Below the dotted line, the lightest right-handed sneutrino corresponds to the real 
component, whereas the imaginary component becomes lighter in regions above that line. 
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sneutrinos are relevant. From the stationary condition 
dV 



dN 



m 2 RR + 2|A JV | 2 |iV| 2 + (2X N (~XH 1 H 2 + kS 2 )* + 2X n A Xn S) ^~ 



= 2X N {-XH 1 H 2 + kS 2 )*N + 2\X N \ 2 \N\ 2 N* 

+ {\2X N S\ 2 + m% + \y N H 2 \ 2 ) N* + 2X n A Xn SN 

= N 

= 0, (2.18) 
the possible minimum or maximum is given by 
~ J 

\\ N )\ =) \2\ N {-\H 1 H 2 +KS ,2 )*+2\ N A XN S\- m 2 RR f f ^ n \ ■ ( 2 - 19 ) 

I 2\x N \ 2 y 11 > U J 

If 2|m|j fl | — iTi 2 Rk < 0, then, a non vanishing VEV is absent and the origin N = 
is the true minimum and vacuum. Of course, this condition is exactly same as the 
condition for the positiveness of the mass squared of the lighter right-handed sneutrino 
in Eq. (2.14). Then, the scalar potential is reduced to the Higgs potential in the 
NMSSM. 

As stated in Introduction, the domain wall problem in this simplest Higgs potential 
can be solved by some modifications [33, 34], which have no impact on the sneutrino 
physics we are interested in. 



2.4 Higgs and neutralino sectors 

The inclusion of the S superfield implies two new Higgs states, one CP-even and one 
CP-odd, which mix with the MSSM Higgs states. We assume that there is no CP- 
violation in the Higgs sector and therefore CP-even and CP-odd states do not mix. 
Thus the corresponding mass matrices, which have exactly the same expressions as in 
the NMSSM, can be written in the respective basis, and the mass eigenstates expressed 
as a linear superposition of the corresponding gauge eigenstates. For CP-even and CP- 
odd Higgses, we have 

Hi = $Hf HuR + H® R + SffO S R , (2.20) 

and 

A° a = P\« H° uI + P 2 o H° dI + Pi„ S t , (2.21) 
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for i — 1, 2, 3, a — 1, 2, and where the subscripts R and I denote the real and imaginary 
components of the corresponding fields, respectively. The minimisation of the scalar 
potential and subsequent calculation of the Higgs masses is carried out with the code 
NMHDECAY 2.0 code [41]. 

The superpartner of the singlet Higgs, the singlino S, has the same quantum num- 
bers as the bino, wino and Higgsinos and therefore mixes with them, giving rise to a 
fifth neutralino state. This is exactly the same situation as in the NMSSM [34, 36, 30] 
and the neutralino mass matrix has the same expression. The mass matrix is diagonal- 
ized as usual by an unitary matrix TV*. The neutralino mass eigenstates are therefore 
a linear superposition of bino, wino, Higgsinos and singlino which we express as 

X° = Nl.B + N?.W 3 + JV|.tf? + jV|.#° + N^S. (2.22) 
2.5 Right-handed sneutrino LSP 

The phenomenology of this construction is largely dependent on which particle plays the 
role of the lightest supersymmetric particle. In the ordinary NMSSM the neutralino is 
the LSP in extensive regions of the parameter space, thereby providing an interesting 
candidate for dark matter, which can reproduce the correct relic abundance and be 
within the sensitivity of some of the future direct detection experiments [30]. 

In this extended model another possibility arises, since the right-handed sneutrino 
might also be the LSP. The mass of the lightest right-handed sneutrino, is given by 
m RR — 2\m RR \ (2.13), and is therefore dependent on the set of new parameters. In 
particular, as we have shown in Fig. 2, below a certain value of \n it is always possible 
to choose a value of the soft mass parameter, m^, for which the right-handed sneutrino 
mass will be small enough so as to guarantee that it becomes the LSP. Being a neutral 
particle with electroweak scale interactions, the right-handed sneutrino would therefore 
constitute a dark matter candidate within the category of WIMPs. We will explore 
this possibility by analyzing the sneutrino relic abundance and direct detection cross 
section in the following Sections. 

3 Thermal relic density 

A crucial feature of this scenario is the existence of direct couplings of the right-handed 
sneutrino to Higgses and neutralinos. These emerge through the term XjySNN in the 
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superpotential (2.2), since as we stated above, the singlet and singlino components of S 
mix with the CP-even Higgs bosons and neutralinos, respectively. The strength of the 
interaction is therefore dependent on the value of A at and A\ N . This provides tree-level 
interactions with ordinary matter which are specified in Appendix A. For adequate 
values of Aat (and A\ N ), these couplings would be of electroweak scale, thereby making 
the right-handed sneutrino a potential WIMP candidate 4 . Still, in order to determine 
whether or not the right-handed sneutrino is a viable WIMP, its relic density has 
to be evaluated and compared with the recent measurement by the WMAP satellite, 
0.1037 < Vlh 2 < 0.1161. [37]. 

The possible annihilation products for right-handed sneutrinos in this construction 
include the following channels 

(i) W + W~, Z Z, and // via s-channel Higgs exchange; 

(ii) H® H®, via s-channel Higgs exchange, t- and -u-channel sneutrino exchange, and a 

scalar quartic coupling; 

(iii) A° a and Hf H~ , via s-channel Higgs exchange, and a scalar quartic coupling; 

(iv) Z A° a and W ± if ^ via s-channel Higgs exchange; 

(v) iViV, via s-channel Higgs exchange and via t- and w-channel neutralinos exchange. 

The processes suppressed by the neutrino Yukawa (such as s-channel sneutrino 
annihilation mediated by the Z boson) have not been included, since these are negligible 
as shown in the previous section. The corresponding expressions for the amplitudes of 
each channel are listed in Appendix B. 

Since all the above annihilations involve s-channel Higgs exchange processes, res- 
onant sneutrino annihilation will take place near the pole, where ~ 2m H g. This 
implies that the partial-wave expansion method is not sufficiently accurate. Thus, we 
integrate it directly, following the same procedure detailed in Ref. [38] for the case of 
the neutralino. In our calculation, we assume the absence of a degenerate particle 
and do not include possible associated co-annihilation effects. Notice that this could in 
principle happen between the real and imaginary component of sneutrinos and between 
the lightest sneutrino and the lightest neutralino. The former situation only arises if 

4 This is contrary to the case of right-handed sneutrinos with the MSSM, which do not present 
(unsuppressed) tree-level couplings to ordinary matter and therefore have to be generated via non- 
thermal processes (e.g., late decays of the NLSP) [15, 16, 17, 18]. 
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m RR in (2.14) is very small, which implies some accidental cancellation in (2.13) and 
only occurs for specific choices if input parameters which we will explicitly avoid. The 
latter case would only be relevant near those regions in which the sneutrino mass is 
close to the NLSP mass. These regions of the parameter space can be easily identified 
and we will comment on this later. 

As we explained above, and is evidenced in the different annihilation channels 
previously detailed, the main feature of this construction is the direct coupling of 
the sneutrino to the Higgs fields. Hence, the sneutrino annihilation cross section is 
extremely dependent on the structure of the Higgs sector. Although this introduces 
a strong dependence of our results on the NMSSM parameter space, there are some 
general features which are easy to identify and understand. 

• The coupling A a? determines the overall scale of the annihilation cross section. 
A larger A a? implies more effective sneutrino annihilation, and in turn a smaller 
relic abundance, and viceversa. Regarding the supersymmetric spectrum, notice 
that Aat only affects the right-handed neutrino and sneutrino masses and does 
not alter the rest of the NMSSM spectrum. Thus, having chosen a set of viable 
NMSSM input parameters, the value of A a? (as well as and A\ N ) can be 
freely varied in order to reproduce the correct sneutrino relic abundance, being 
only constrained by perturbativity (we will impose A^ < 1), the occurrence of 
tachyons in the sneutrino sector or the upper bound on the sneutrino mass in 
order for it to be the LSP. 

• Annihilation into H® H® and A Q a A® will turn out to be among the most effective 
channels, as we will later see. Whether these are kinematically allowed or not 
depends on the Higgs masses. Interestingly, within the framework of the NMSSM, 
very light CP-even and CP-odd Higgses are possible (as long as they have a 
significant singlet component), making these channels available for a wide range 
of sneutrino masses. In the next Section we will show how this can, for instance, 
make it possible for light sneutrinos to reproduce the correct relic abundance. 

• Another relevant contribution to the total annihilation cross section is due to 
the annihilation into a pair of right-handed neutrinos, N. This channel is kine- 
matically allowed when > Mjv- As we saw in Fig. 1, this is allowed above 
a certain value of the soft sneutrino mass, m^, which depends on the specific 
values of Aat and A\ N and which increases for larger \ N . 
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• Finally, for all the possible annihilation products there is always a contribution 
coming from s-channel CP-even Higgs exchange. This implies that all of them are 
subject to a resonant effect when 2m ^ ~ f^u ) f° r * = 1,2,3. The occurrence of 
resonant sneutrino annihilation gives rise to a characteristic decrease of the relic 
abundance at the corresponding values of the sneutrino mass. Furthermore, given 
the possibility of light scalar Higgses in the NMSSM, this resonant annihilation 
can be present even for light sneutrinos. 

Having understood the basic properties of right-handed sneutrino annihilation in 
this framework, we will now proceed to explore the parameter space and provide specific 
examples which show how the correct thermal relic abundance can be obtained. 

3.1 Numerical examples 

The input parameters of this scenario are, on the one hand, the usual NMSSM param- 
eters, which we define at low-energy, 

A, k, tan/3, n, A\, A K . (3.23) 

The soft supersymmetry-breaking terms, namely gaugino masses, M 1j2j3 , flavour inde- 
pendent scalar masses, ^q,l,i/,l>,.e, and trilinear parameters, Au^^-, are also taken as 
free parameters and specified at low scale. More specifically, we have assumed they 
gaugino masses mimic, at low-energy the values obtained from a hypothetical GUT 
unification. Thus we set M 1 : M 2 : M 3 = 1:2:6. The present lower bound on a 
possible super symmetric contribution to the muon anomalous magnetic moment, sets 
stringent upper constraints on the mass of sleptons. As it was shown in [30], for small 
values of tan (5 the experimental result can be reproduced with the choice m^E = 150 
GeV and A E = —2500 GeV and a small value for the gaugino mass 5 , Mi < 160 GeV. 
Larger values of tan/3 are, however less constrained [39]. In order to satisfy the exper- 
imental constraint on the branching ratio of the b — > 57 rare decay a careful choice has 
to be made of parameters in (3.23). More details on the conditions under which these 
bounds can be fulfilled can be found in [30, 40]. 

5 It should be stressed that if one does not wish to impose the bound on the muon anomalous 
magnetic moment (motivated e.g. by tau data, which lead to a better agreement with the SM result), 
large slepton masses, equal to squark masses can be chose. This would, however, have no consequence 
on the calculation of neither the relic abundance of right-handed sneutrinos, nor in their detection 
cross section. 
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The analysis of the low-energy NMSSM phenomenology has been performed with 
the NMHDECAY 2.0 code [41], which minimises the scalar potential, dismissing the pres- 
ence of tachyons and/or false minima, and computes the Higgs boson masses including 
1- and 2-loop radiative corrections, as well as the rest of the supersymmetric masses. 
Based on this code, we have built a set of routines which numerically calculate the 
right-handed sneutrino spectrum and relic density as described in the previous section. 

More specifically, we impose the experimental bound on the branching ratio of the 
rare b — > S7 decay, 2.85 x 10~ 4 < BR(b — > sj) < 4.25 x 10~ 4 at 2a level, obtained 
from the experimental world average reported by the Heavy Flavour Averaging Group 
[44], and the theoretical calculation in the Standard Model [45], with errors combined 
in quadrature. We also take into account the upper constraint on the Bs — > fi + fi~ 
branching ratio obtained by CDF, BR(B S -> < 5.8 x 1(T 8 at 95% c.l. [42] 

(which improves the previous one from DO [43]). 

Regarding the muon anomalous magnetic moment, a constraint on the supersym- 
metric contribution to this observable, a^ USY , can be extracted by comparing the ex- 
perimental result [46], with the most recent theoretical evaluations of the Standard 
Model contributions [47, 48, 49]. When e + e~ data are used the experimental excess 
in a M = (g^ — 2)/2 would constrain a possible supersymmetric contribution to be 
a^ USY = (27.6 ± 8) x 10~ 10 , where theoretical and experimental errors have been com- 
bined in quadrature. However, when tau data are used, a smaller discrepancy with the 
experimental measurement is found. Due to this reason, in our analysis we will not 
impose this constraint, but only indicate whether the points are compatible with it at 
the 2a level, for which we assume the range 11.6 x 10~ 10 < a^ USY < 43.6 x 10~ 10 . 

The inclusion of the new superfield N and the corresponding terms in the su- 
perpotential and Lagrangian leaves three new parameters to be fixed. Following the 
discussion in the previous section, these can be chosen as 

Ajv, m^, A Xn . (3.24) 

These, together with the NMSSM parameters (3.23) fully specify the model. Our task 
is now to determine whether the correct sneutrino relic density can be obtained with 
reasonable choices of the above set of parameters. 

Let us now illustrate with some examples the theoretical predictions for the relic 
abundance of right-handed sneutrinos and the various possibilities which allow to re- 
produce the correct relic density. The sets of parameters for each example are detailed 
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A) 


Bl) 


B2) 


C) 


tan/3 


5 


5 


5 


3.5 


A x 


550 GeV 


400 GeV 


400 GeV 


480 GeV 


A K 


-200 GeV 


GeV 


GeV 


-50 GeV 


n 


130 GeV 


200 GeV 


200 GeV 


230 GeV 


(A, «) 


(0.2, 0.1) 


(0.1, 0.05) 


(0.3, 0.2) 


(0.5, 0.3) 


Mi 


200 GeV 


150 GeV 


150 GeV 


300 GeV 




250 GeV 


250 GeV 


250 GeV 


250 GeV 




1000 GeV 


1000 GeV 


1000 GeV 


1000 GeV 


^£ 


-2500 GeV 


-2500 GeV 


-2500 GeV 


-2500 GeV 


Au,d 


1500 GeV 


1500 GeV 


1500 GeV 


1000 GeV 


m H o 


62.4 GeV 


117.7 GeV 


116.2 GeV 


115.1 GeV 


m H o 


119.4 GeV 


200.8 GeV 


267.7 GeV 


263.5 GeV 


m H o 


634.1 GeV 


706.4 GeV 


730.6 GeV 


725.4 GeV 


m A o 


199.6 GeV 


14.2 GeV 


47.6 GeV 


169.7 GeV 


m A o 


632.5 GeV 


705.1 GeV 


727.9 GeV 


719.5 GeV 




95.8 GeV 


126.3 GeV 


124.3 GeV 


190.1 GeV 


BR(b -> S7) 


4.15 x 10~ 4 


4.02 x 10~ 4 


4.02 x 10~ 4 


4.03 x 10~ 4 


„SUSY 


4.08 x 10~ 10 


5.58 x 10- 10 


8.05 x 10~ 10 


2.96 x 10- 10 



Table 1: Set of inputs corresponding to the examples used in the analysis, together with 
the predicted values of BR(b — > S7) and a^ USY and the relevant part of the resulting 
supersymmetric spectrum. 
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in Table 1, together with the predicted values 6 of BR(6 — > 57) and a^ USY . 
3.1.1 Predominant annihilation into scalar Higgses 

As already mentioned, sneutrino annihilations into Higgs bosons are generally among 
the dominant channels, when they are kinematically allowed. This is also the case 
in other models for thermal right-handed sneutrino [23]. However, in our case, the 
flexibility and interesting properties of the Higgs sector of the NMSSM make this 
possibility much very versatile. 

As an specific example we will choose the NMSSM input parameters as in case A) of 
Table 1. This point is very characteristic of the NMSSM since light CP-even Higgses are 
possible. In this case the lightest Higgs has a mass m H o ps 62 GeV, which is consistent 
with the present LEP constraints due to its large singlet component ((5^o) 2 = 0.997). 
The second lightest Higgs is MSSM-like and has a mass of m H o ps 120 GeV. The lightest 
neutralino is a mixed bino-Higgsino-singlino state ((A^-J 2 ps 0.1, (iV-J 2 + (A"|J 2 ps 0.7 
and (iV|J 2 p^ 0.2). Its mass is m^o ps 96 GeV, and this sets the upper bound for the 
sneutrino mass if the latter is to be the LSP. Regarding the right-handed sneutrino 
sector, its mass is a function of the parameters A at, and A\ N , as explained in the 
previous section. As a first example we have fixed A\ N = —250 GeV and performed a 
scan in the allowed range for the A at, parameters. 

The resulting sneutrino annihilation cross section is shown in Fig. 4, together with 
the contribution for each individual channel, as a function of the sneutrino mass for 
A at = 0.06 (the variation in the sneutrino mass is achieved with a scan in m^). A first 
thing to notice is that the various annihilation channels display a resonant enhancement 
due to the two lightest CP-even Higgses when the sneutrino mass is approximately one 
half of the Higgs mass for the two lightest scalar Higgses, i.e., « 30 GeV and 
60 GeV. The resonance due to the third scalar Higgs state is not observed since it 
would occur for large masses (m^ ps 317 GeV), for which the right-handed sneutrino 
is no longer the LSP. 

In this example, annihilation into a pair of lightest CP-even Higgses dominates the 
total annihilation cross section once this channel becomes kinematically allowed (for 
irifj^ > m H o ~ 62 GeV). Lighter sneutrinos can only annihilate into a pair fermion- 

6 Some of these examples are inspired in the scenarios analysed in [30] (e.g., cases Bl) and B2) are 
similar to the case studied in Fig. 7 there), and we refer the reader to that work for a more detailed 
study of the NMSSM parameter space and associated experimental constraints. 
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Figure 4: Left) Total sneutrino annihilation cross-section as a function of the sneutrino mass 
(grey thick solid line) for example A) of Table, 1 with (X,k) = (0.2,0.1), and A^v = 0.06. 
The contributions from the different annihilation channels are indicated as follows, dashed 
red line for //, dot-dashed green lines for ZZ and W + W~ , solid blue line for H®Hj, and 
solid black line for iViV. Vertical dashed lines indicate the location of the resonances in the 
s-channels when 2m^ ~ m H o, whereas vertical dotted lines indicate when the different 
annihilation channels become kinematically accessible. The vertical solid line indicates the 
point at which the neutralino becomes the LSP. Right) Theoretical predictions for the 
sneutrino relic density as a function of the sneutrino mass for the same example, but for a 
scan with Ajv £ [0.05,0.1] and and 6 [0,150] GeV. The solid line indicates the result 
for X N = 0.06. 



antifermion, being bb the leading contribution. As observed in the plot, for this choice 
of parameters this channel is less effective and only provides a sufficiently large an- 
nihilation cross section close to the resonances of the scalar Higgses 7 . Notice also 
that although annihilation into a pair of right-handed neutrinos is also possible for 
sneutrinos heavier than 80 GeV, this channel only becomes comparable to the Higgs 
contribution for sneutrino masses of order 150 GeV, for which the sneutrino is no longer 
the LSP. 

The resulting theoretical predictions for the relic density are depicted on the right 

7 We remind the reader that for a WIMP of mass 100 GeV the necessary annihilation cross-section 
in order to reproduce the WMAP result is around a ann ~ 0.1 pb. 
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Figure 5: Effect of the relic density constraint on the (\ N ,m N ) P' ane f° r case A) in 
Table 1 with A\ N = —250 GeV. Black dots represent the regions where the sneutrino 
relic abundance is in agreement with the WMAP constraint, whereas light(dark) grey dots 
represent those where the sneutrino relic abundance is smaller(larger). The curved solid 
lines indicate the trajectories with a fixed value of the sneutrino soft mass m N . 

hand-side of Fig. 4 for a scan of A^ in the range Xn £ [0.05,0.1] and a scan in the 
allowed values for m N . The solid line indicates the result for Xn = 0.06. We can see 
how light sneutrinos with masses as small as approximately 25 GeV can be obtained 
with the correct relic abundance due to the resonant effects in the Higgs diagrams. 
Once annihilation into scalar Higgses is allowed, the correct relic abundance can be 
obtained in the whole range 8 m Ni ^ 60 — 90 GeV with \ N 0.06. 

For a better understanding of the effect of the relic density constraint on the sneu- 
trino parameters, we represent in Fig. 5 the (Xn, m N 1 ) plane, indicating with black dots 
the regions where the sneutrino relic abundance is in agreement with the WMAP result. 
Dark grey dots represent the areas in which the sneutrino relic abundance exceeds the 
WMAP constraint and light grey dots are those in which the sneutrino relic density is 
smaller than the WMAP bound. The contours for a fixed sneutrino soft mass, m N , are 
indicated by means of solid lines. As already explained, when annihilation into a pair 
of Higgs bosons is allowed, a value of Xn ~ 0.06 is favoured in this example. As we 

8 Since co-annihilation effects have not been included, our result for the relic abundance is not exact 
when mf Ji w m^o . 
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m N1 (GeV) m N1 (GeV) 

Figure 6: Left) The same as in Fig. 7 (case A) in Table 1) but with A\ N = —500 GeV 
and \ N = 0.25. Right) Theoretical predictions for the sneutrino relic density as a function 
of the sneutrino mass for the same example, but for a scan with Xjy e [0.05,0.35] and 
and mfj e [0,250] GeV. The solid lines indicate, from top to bottom, the result for 
A^v = 0.11, 0.25. 

see, this corresponds to values of the sneutrino soft mass in the range 100 — 150 GeV. 
Otherwise the relic density is too large except along the resonances (the fine-tuned 
near horizontal strips), for which the value of A^ can be as large as 0.15. 

3.1.2 Predominant annihilation into bb and very light sneutrinos 

Obtaining the correct relic density with only annihilation into bb would require a larger 
value of Xn- However, as we can see from Fig. 5, if we restrict ourselves to positive values 
of mjy then for this choice of parameters there is an upper bound A at < 0.15 — 0.18. 
Remember however from the previous section that an increase in \A\ N \ allows larger 
values of A at- With this in mind, we have reanalysed the previous example, case A) in 
Table, 1, but now using A\ N = —500 GeV, for which values of A^v up to 0.35 can be 
used. This considerably increases the annihilation cross section. 

For this example, Fig. 6 represents the annihilation cross section for A^ = 0.25 and 
a scan in the allowed values of m^. The bb contribution is now significantly larger, 
reaching values of order of 1 pb without the need of resonant effects. Similarly, the Higgs 
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contribution also increases considerably, thereby leading to a too large annihilation 
cross section (notice the different scale on the plot). On the right-hand side of Fig. 6 
the resulting sneutrino relic abundance is plotted for a scan with \jy £ [0.05,0.35] 
and and G [0,250] GeV, where the solid lines indicate, from top to bottom, the 
particular cases \n = 0.11 (which yields similar results as the example studied in 
Section 3.1.1) and 0.25. Interestingly, as we can see, for such large values of \ N light 
sneutrinos with masses > 5 GeV can reproduce the correct relic abundance. 

3.1.3 Predominant annihilation into pseudoscalar Higgses 

Another interesting possibility is sneutrino annihilation into a pair of pseudoscalar 
Higgses. Although the pseudoscalar is rather heavy in the MSSM (and experimentally 
constrained to be above 93.4 GeV [62]), in the NMSSM (where an extra CP-odd state 
is present) it can be much lighter in some regions of the parameter space and experi- 
mentally allowed if it has a large singlet component [63, 64, 65]. This therefore implies 
that this channel can be kinematically allowed for lighter sneutrinos. 

In order to illustrate this possibility, we have taken the NMSSM input parameters 
as indicated in example Bl) of Table 1. This point in the NMSSM parameter space is a 
good example of how light pseudoscalars might be phenomenologically viable. In this 
case the lightest pseudoscalar has a mass m A o = 14.2 GeV and it is singlet-like. The 
lightest scalar Higgs is MSSM-like with a mass m H o « 116 GeV, whereas the second 
lightest Higgs, with m H o 200 GeV is mostly singlet. The lightest neutralino is a 
mixed bino-Higgsino state ((NiJ 2 « 0.66 and (-/VfJ 2 + (N^) 2 « 0.25) with a mass 
m^o = 126 GeV. It would be the lightest supersymmetric particle in the NMSSM for 
this choice of parameters and therefore sets once more the upper limit for allowed values 
of the right-handed sneutrino mass. Regarding the right-handed sneutrino sector, we 
have fixed A\ N = —500 GeV and performed a scan in the A at, nifj parameters. This 
is precisely the same example that was displayed on Fig. 2. As one can extract from 
there, a maximal value of Ajv ~ 0.11 can be used, with the soft sneutrino mass ranging 
from = to approximately = 235 GeV. 

The predicted sneutrino annihilation cross section is represented on the left hand- 
side of Fig. 7 as a function of the sneutrino mass, together with the contributions from 
each individual annihilation channel. As we can see, for this example the annihilation 
into a pair of CP-odd Higgses largely dominates over the other contributions. This is 
particularly true for light sneutrinos since, due to the smallness of m A o, this channel 
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Figure 7: Left) The same as in Fig. 4 but for example Bl) of Table, 1 with A = 0.1, k = 0.05, 
and \ N = 0.07. Here the dashed blue line denotes the contribution of the A® A® channel and 
the dotted blue line (only visible in the Higgs resonances) the contribution of A\Z . Right) 
Theoretical predictions for the sneutrino relic density as a function of the sneutrino mass 
for the same example, but for a scan with G [0.05,0.1] and and E [0,250] GeV. 
The solid line indicates the result for Aat = 0.07. 

is kinematically available for > 14 GeV. Note that for lighter sneutrinos the only 
possibility is annihilation into a pair of //, which is less effective. The abrupt enhance- 
ments at the resonances of the Higgs-exchanging s-channels when 2m ^ ~ m H o,m H o 
are also observed for every channel at ~ 58 GeV and 100 GeV. 

The theoretical predictions for the resulting relic density are depicted on the right 
hand-side of Fig. 7 for the range Ajv £ [0.05,0.11]. As we see there, the correct relic 
density can be obtained for a wide range of sneutrino masses. In particular, sneutrinos 
as light asms ~ 15 GeV are now possible due to their very effective annihilation into 
a pair of CP-odd Higgses. This is particularly interesting, since these light sneutrinos 
have different properties from those obtained in case A). In case A) they could only 
annihilate into bb and thus the coupling A at had to be of order 0.25. However, in the 
present case the coupling can be much smaller and Aat ~ 0.07 suffices to obtain the 
correct relic abundance. This will have very interesting properties for their detection, 
as we will see in Section 4. 
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Figure 8: The same as in Fig. 5 but for case Bl) of Table 1 with A Xn = —500 GeV. 

The effect of the relic density constraint on the (Ajy, ) plane is shown in Fig. 8, 
together with the contours for a fixed sneutrino soft mass, m^. As we observe, values 
of the sneutrino coupling constant to the Higgs of order Ajv ~ 0.06 — 0.1 are sufficient 
to have viable sneutrino dark matter, whereas values of the soft sneutrino mass around 
mfj pa 100 — 200 GeV seem to be preferred. Once more, along the Higgs resonances 
sneutrino annihilation is more effective and smaller values of A^ are needed to obtain 
the correct abundance. Remarkably, the region corresponding to viable very light 
sneutrinos (with masses in the 15 — 40 GeV range) can be obtained without requiring 
a fine-tuning in the A at or parameters 9 . 

The sneutrino properties are extremely sensitive to the NMSSM Higgs sector. One 
should therefore expect that the above analysis and preferred values for the sneutrino 
parameters could vary for a different set of NMSSM inputs. In particular, given a 
different point of the (A, k) plane and due to the changes in the (scalar and pseu- 
doscalar) Higgs masses and mixings, the values of the sneutrino mass at which the 
various channels open and the conditions for resonant annihilation would differ. 

This is indeed the case. For instance, for the same choice of parameters as in 

9 Remember, however that there is a certain fine tuning in the NMSSM parameters in order to 
obtain very light viable pseudoscalars. Still, very light right-handed sneutrinos in this construction 
are much more natural than very light neutralinos in the NMSSM, for which the extra condition of a 
resonant annihilation with the very light CP-even Higgs, 2m^o ss m A a is necessary [64]. 
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Figure 9: Left) The same as in Fig. 7 but for case B2) of Table 1 with A^ = 0.14. Right) 
Corresponding relic density as a function of the sneutrino mass for A^ G [0.05,0.35]. The 
solid line indicates the result for Aat = 0.14. 



example the previous example, let us now choose (A, k) = (0.3,0.2), corresponding to 
case B2) in Table 1. This implies an increase of the pseudoscalar mass (and a mild 
reduction of its singlet component) for which m A o « 48 GeV, as well as an increase 
in the mass of the second lightest scalar Higgs, which remains being singlet-like (the 
lightest scalar Higgs is still MSSM-like and has a mass m H o ~ 116 GeV). The rest of 
the spectrum is detailed in Table 1. 

As a consequence of these changes, annihilation into pseudoscalars is no longer 
possible for very light sneutrinos, which can now only do it into bb and their annihilation 
cross-section becomes much smaller (for the same value of Ajv) than in example A). For 
sneutrinos with masses above 50 GeV the pseudoscalar channel is still the dominant 
one, with a modest contribution from annihilation into a pair of gauge bosons when 
these are kinematically allowed. This can be observed on the left-hand side of Fig. 9, 
where the corresponding predictions for the sneutrino annihilation cross section are 
shown for Aa? = 0.14. Notice also that due to the change in the CP-even Higgs masses, 
the position of the resonances change, in particular, the resonance due to the second 
heaviest Higgs is not present (it would occur for ~ 200 GeV, a mass for which 
the sneutrino is no longer the LSP). 
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Figure 10: The same as in Fig. 5 but for case B2) of Table 1 with A\ N = —500 GeV. 

As illustrated on the right-hand side of Fig. 9, the correct relic abundance can still 
be obtained for a given range of sneutrino masses. As in the previous examples, we 
perform a scan in the allowed values of Ajy and for A\ N = —500 GeV. As explained 
in the previous section, due to the increase in A\ N , larger values of Aat are viable. 

The effect of the relic density constraint on the (\ N ,m I ^ i ) plane is displayed in 
Fig. 10. The allowed areas for sneutrinos heavier than 60 GeV correspond to Xn ~ 
0.12 — 0.14 with a sneutrino soft mass of order ps 200 — 225 GeV, whereas lighter 
sneutrinos are only viable in the thin region corresponding to resonant annihilation 
with the lightest scalar Higgs. 

As we see from examples Bl) and B2), through adequate variations in the sneutrino 
parameters (3.24), it is possible to find different points in the (A, k) plane where the 
sneutrino has the correct relic density. Variations on the A and k parameters have 
a large impact on the neutralino and Higgs sector and it is therefore not surprising 
that as a consequence the resulting sneutrino phenomenology is significantly altered. 
However, the flexibility in the sneutrino properties (especially the freedom in choosing 
an adequate value of its couplings through A^) makes it possible to obtain viable 
sneutrino dark matter for a wide choice of NMSSM parameters. 

In order to illustrate this we have combined a scan in the (A, k) plane, in the range 
A, k — 0.07 — 0.7 with a scan in the sneutrino parameters taking = — 250 GeV 
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Figure 11: Left) The same as in Fig. 5 but for a scan in (A, k) as described in the text. Empty 
circles correspond to points with the correct relic abundance due to resonant annihilation 
effects. Right) Regions with the correct sneutrino relic abundance (black dots) in the (A, k) 
plane. The oblique ruled area is excluded due to the occurrence of Landau Poles or because 
of experimental constraints in the Higgs sector, whereas the vertical ruled area is ruled out 
due to tachyons in the Higgs sector. 



and Ajv = 0.05 — 0.3. We have fixed A\ N = —500 GeV and we have taken otherwise the 
same NMSSM parameters as in example Bl) and B2) of Table 1. On the left-hand side 
of Fig. 11 we represent the resulting sneutrino mass as a function of the coupling A^. 
Black dots represent those with the correct relic abundance and light points correspond 
to those where the sneutrino relic density is smaller than the WMAP results (points 
with a relic density exceeding the WMAP constraint are not shown). As we can see, 
the correct relic abundance can be obtained within a range of Ajv ~ 0.06 — 0.15 and 
for a wide range of sneutrino masses. On the right-hand side of Fig. 11 the regions 
with the correct relic abundance are displayed on the (A, k) plane. We clearly see that 
the sneutrino can be a viable candidate (for an adequate choice of the A^ and 
parameters) in virtually any point in this plane. Notice that the absence of points in 
the left-hand side of the plot is due to the fact that the imaginary component of the 
sneutrino becomes the lightest one, as explained in Section 2. This is in principle no 
problem and in fact the imaginary sneutrino would have almost the same properties 
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Figure 12: Left) The same as in Fig. 7 but for case C) of Table 1 with X N = 0.15 and 
A\ N = —250 GeV. Right) Corresponding relic density as a function of the sneutrino mass 
for for Aat G [0.1,0.3]. The solid line indicates the result for Ajv = 0.15. 

as the real component 10 , however, as we mentioned earlier, for simplicity we have not 
considered such a possibility here. 



3.1.4 Predominant annihilation into right-handed neutrinos or gauge bosons 

Let us finally illustrate with another example the relevance of other annihilation chan- 
nels, in particular the annihilation into a pair of right-handed neutrinos. Our set of 
input parameters is chosen according to example C) in Table 1, where we also specify 
the resulting spectrum. 

We fix A\ N = —250 GeV and plot the resulting annihilation cross section for each 
individual channel on the left-hand side of Fig. 12. Once more, the features of the 
Higgs sector clearly affect the calculation of the relic density for the sneutrino. In this 
example the resonant annihilation through the exchange of the two lightest Higgses is 
visible at « 58 GeV and 130 GeV. The lightest Higgs in this example is doublet- 
like, whereas the intermediate Higgs, H2, and the lightest pseudoscalar are almost 

10 The only difference arises in the expression for its annihilation into right-handed neutrinos, as 
explained in Appendix B. 
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Figure 13: The same as in Fig. 5 but for case C) of Table 1 with A\ N = —250 GeV. 

pure singlets. Sneutrino annihilation into bb is the only possibility for sneutrino masses 
below Mw, but this channel is only effective when resonant annihilation through the 
lightest Higgs takes place. Annihilation into ZZ or W + W~ is the dominant channel 
below the Higgs production threshold. Interestingly, after the second Higgs resonance, 
annihilation into a pair of right-handed neutrinos becomes the leading contribution for 
m N 1 ~ 1^0 GeV, although the contributions of the H®H® and A® A® channels are not 
negligible. 

The resulting relic density is shown on the right-hand side of Fig. 12 for a scan with 
An G [0.1,0.3] and G [0,200] GeV. The solid line indicates the result for A^ = 
0.15. We find that in this example the correct relic density can be obtained through 
annihilation into ZZ or W + W~ (for sneutrinos with a mass around 100 GeV and 
\n ~ 0.2) or annihilation into NN and H®H® for sneutrinos above > 140 GeV, 
in which case A^ w 0.1 is necessary. 

The corresponding (Ajy, rrijj ) plane is shown in Fig. 13, indicating the regions con- 
sistent with the WMAP constraint on the relic density. Consistently with the discussion 
above, in the region above the dot-dashed line (which corresponds to the area in which 
the sneutrino mass is larger than the right-handed neutrino mass, m^> Mat) there 
is a wide region with A at ~ 0.1 for which the WMAP constraint is satisfied. This also 
corresponds to values of the soft sneutrino masses in the range = 150 — 200 GeV. 
Notice also a second region for larger values of A^ of order 0.15 — 0.18. Below the 
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dot-dashed line, where < the WMAP region corresponds to values of the 
sneutrino-Higgs coupling in the range Ajv = 0.15 — 0.24, and smaller values of the 
sneutrino soft mass, = 100 — 150 GeV. Light sneutrinos in this example have a too 
large relic abundance and it is only along the regions with resonant annihilation that 
the WMAP constraint is satisfied. 

3.2 Overview 

Throughout this section we have seen how the right-handed sneutrino can reproduce 
the correct dark matter abundance under different conditions and for various points in 
the NMSSM parameter space. In this sense, the flexibility to choose adequate values of 
the sneutrino parameters A at, and A\ N was very useful. Notice in any case that the 
preferred values of these parameters are quite natural. The sneutrino-Higgs coupling 
\n can be chosen in the range \ N = 0.05 — 0.4 and the sneutrino soft mass could be 
of order of 100 — 250 GeV, being of the same order than the slepton soft masses in our 
examples. As a result, sneutrinos in the mass range 5 — 200 GeV are found 

that can reproduce the WMAP constraint on the relic abundance. 

It should also be noted that obtaining heavier right-handed sneutrinos is in principle 
possible, especially since more channels would be kinematically allowed and the relic 
abundance could be obtained for moderate values of A at. Notice that, in order to do so, 
the mass of the neutralino (and that of the rest of the spectrum) also must be increased 
so that the right-handed sneutrino is the LSP. This can be done through an increase 
in the gaugino and scalar mass parameters as well as in the \i term. However, this 
entails a further reduction of the supersymmetric contribution to the muon anomalous 
magnetic moment, thereby increasing the tension with its experimental bound. 

4 Direct detection 

The direct detection of sneutrinos could take place through their elastic scattering with 
nuclei inside a dark matter detector. The interaction occurs in the non-relativistic 
regime (given that the velocity of sneutrinos in the dark matter halo is very small) 
and therefore one can easily find a description in terms of an effective Lagrangian. 
In our case, there is only one diagram contributing (at tree level) to this process, 
namely, the t-channel exchange of neutral Higgses shown in Fig. 14 (the exchange of 
a Z boson is largely suppressed by the neutrino Yukawa squared and therefore is 
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Figure 14: Diagram describing the elastic interaction of sneutrinos with quarks. 

completely negligible). The effective Lagrangian describing the four-field interaction 
only contains a scalar coupling which reads 

C eff D a qi NNq iqi (4.25) 

with 

lit Ho 



where C H o^ is defined in Appendix A, Y qi is the corresponding quark Yukawa coupling 
and i labels up-type quarks (i = 1) and down-type quarks (i = 2). Notice that the 
effective Lagrangian contains no axial-vector coupling since the sneutrino is a scalar 
field, therefore implying a vanishing contribution to the spin-dependent cross section. 
The total spin-independent sneutrino-proton scattering cross section yields 

= f P , 2 f P , (4-27) 

1Vp 7r {m p + m^Y p 

where m p is the proton mass and 



m p ^ Tqt m qi 27 ' * " <— ' m q% 

' qi=u,a,s qi=c,b,t ** 



The quantities and /£ G are the hadronic matrix elements which parametrize the 
quark content of the proton. They are subject to considerable uncertainties [66, 67, 68] 
which induce a significant correction to the theoretical predictions for cr^ . In our 
analysis we will consider the most recent values for these quantities, as explained in 
[69]. 

It is obvious from the previous formulae that the sneutrino detection cross section is 
extremely dependent on the features of the Higgs sector of the model. In particular, 
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becomes larger when the sneutrino-sneutrino-Higgs coupling increases (which according 
to its expression in Appendix A can be achieved by enhancing or with large values 
for IAxjJ). Moreover, from Eq. (4.26) we also see that larger values a^f can be obtained 
in those regions of the parameter space where the mass of the lightest Higgs becomes 
smaller. 

We will now calculate the theoretical predictions for the sneutrino scattering cross 
section for the various cases studied in Section 3. Let us begin by addressing case A) in 
Table 1, for which we have seen in Section 3.1.1 that the H®H® channel is the dominant 
contribution to the annihilation cross section for moderate values of Xn- Choosing 
A\ N = —250 GeV and performing the same scan in Xn and m N as in Fig. 4 we have 
calculated the resulting sneutrino scattering cross section off quarks in those cases 
where the sneutrino abundance is consistent with the WMAP constraint or smaller (in 
which case the sneutrino could be a subdominant dark matter component). 

If the sneutrino only contributes to a fraction of the total dark matter density, one 
should expect that it is present in the dark matter halo in the same proportion as in 
the Universe, contributing only to a fraction of the local dark matter density. This 
implies, of course, a reduction of the detection rate in direct detection experiments. In 
order to take this effect into account, it is customary to define the sneutrino fractional 
density £ = min[l, fl N h 2 /0.1037] [50] and plot the modulated elastic scattering cross 
section, fcrf? . 

' ' Np 

We follow this approach and show on the left-hand side of Fig. 15 the theoretical 
predictions for the spin-independent contribution to £cr^ as a function of the sneutrino 
mass for the example studied in Section 3.1.1. This is, case A) with A\ N = —250 GeV, 
Xn £ [0.05,0.1] and m N e [0,250] GeV. Black dots correspond to points with a relic 
density consistent with the WMAP results, whereas grey dots stand for those with 
Vt N h 2 < 0.1037. 

For light sneutrinos, for which annihilation into Higgses is not kinematically al- 
lowed, the value of Xn in order to have the correct relic density is large (see Fig. 5). 
Consequently, the resulting scattering cross section turns out to be sizable and easily 
exceeds the sensitivities of present detectors. On the other hand, when annihilation 
into CP-even Higgses is possible, the necessary value of Xn is smaller, of order 0.1, 
resulting in £,<J^ p < 10~ 8 pb. These results are not yet excluded by present searches, 
but interestingly, could be within the reach of future experiments such as SuperCDMS. 
Notice also that along the Higgs resonances the sneutrino scattering cross section de- 
creases considerably since the necessary value of the Xn coupling to obtain the correct 
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Figure 15: Theoretical predictions for £er^ , as a function of the sneutrino mass for example 
A) of Table 1. All the points represented fulfil all the experimental constraints. Black dots 
correspond to those which reproduce the WMAP result for the dark matter relic density 
whereas grey ones represent those with Vl^h 2 < 0.1037. The sensitivities of present and 
projected experiments are represented by means of solid and dashed lines, respectively, in 
the case of an isothermal spherical halo. The large (small) area bounded by dotted lines is 
consistent with the interpretation of DAMA experiment in terms of a WIMP. 

relic density becomes much smaller. 

On the right-hand side of Fig. 15 we show the results for case A) but for A\ N = 
—500 GeV, which corresponds to the case analysed in Section 3.1.2. We remind the 
reader that in this example a larger A at was allowed and we perform a scan in A at G 
[0.05, 0.35] with 6 [0, 250] GeV. This made it possible to obtain very light sneutri- 
nos with the correct relic abundance. The large value of A^ implies that the detection 
cross section for these light sneutrinos is also very large. This is clearly evidenced in 
the plot, where we find that sneutrinos lighter than 20 GeV have a cross section which 
can be as large as £<r^ ~ 10 -4 pb. 

This situation is very similar to what happens with very light neutralinos in the 
MSSM. Very light neutralinos with masses m^o > 7 GeV can be obtained [51, 52] if 
the GUT relation for gaugino masses is abandoned. Although it was argued that these 
neutralinos could account for the DAM A/LIBRA annual modulation signal without 
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Figure 16: The same as in Fig. 15 but for a scan in the (A, k) plane. Empty circles 
correspond to points with the correct relic abundance due to resonant annihilation effects. 

contradicting the null results from CDMS and XENON10 [53], this interpretation is now 
more constrained by the recent results published by the MAJORANA collaboration 
using results from low-threshold experiments [54]. 

Finally, in order to explore more general variations in the NMSSM parameters, we 
have scanned in the allowed (A, k) plane, keeping A Xn = —500 GeV and the same scan 
in the \n and parameters as in the previous case. As we see from the resulting 
cross section, depicted in Fig. 16 where we see again that the cross section can be rather 
large (especially in those cases with light CP-even Higgs bosons). The region with very 
light sneutrinos, for with the predicted is sizable, could be within the reach of 
low-threshold experiments. 

Next we address the examples Bl) and B2) studied in Section 3.1.3, and for which 
annihilation into pseudoscalar Higgses was predominant. The results for the sneutrino 
scattering cross section are illustrated in Fig. 17. 

We find a very different behaviour from the previous examples. First, since light 
sneutrinos are in this case obtained thanks to their annihilation into very light pseu- 
doscalars, the value of \n at which the WMAP result is reproduced is considerably 
smaller. Second, the mass of the CP-even Higgs bosons is larger in this example than 
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Figure 17: The same as in Fig. 15 but for examples Bl) and B2) of Table 1. 

in the previous one. As a consequence, these very light sneutrinos have a smaller scat- 
tering cross section and can evade the experimental constraints from XENON10 and 
CDMS. 

For example, on the left-hand side of Fig. 17 (which corresponds to the case analysed 
in Fig. 7) sneutrinos as light as ~ 15 GeV are obtained with the correct relic 
abundance and with a scattering cross section of order C,a^ ~ 2 x 1CT 9 pb. This is a 
very interesting situation, since as we will explain later in more detail, it is a completely 
different behaviour than the one observed for neutralinos in the MSSM. 

On the right hand-side of Fig. 17 we show the results corresponding to the case 
analysed in Fig. 9). In this example the CP-even Higgs is heavier than in the previous 
examples and as a consequence we observe a decrease in the scattering cross section. 
In this case ~ 1CT 9 pb and Ton-scale detectors, such as the projected XENON1T 
would be necessary to explore them. 

For completeness we have also performed a scan in the (A, k) plane, the results 
of which are shown on Fig. 18. Depending on the point of the (A, k) plane, the mass 
and coupling of the lightest pseudoscalar varies and the predictions for sneutrino dark 
matter are seriously affected. As a consequence, the predicted cross section for very 
light sneutrinos spans several orders of magnitude, ranging from 1CT 11 pb to 1CT 6 pb. 
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Figure 18: The same as in Fig. 17 but for a scan in the (A, k) plane. 

Interestingly very light sneutrinos populate an area of the a^ — m^ plane which is for- 
bidden to neutralinos in the MSSM. However, the predictions for very light neutralinos 
in the NMSSM are very similar to those obtained here for sneutrinos [54]. 

Finally let us address an example which we already studied in Section 3. 1.4 and for 
which annihilation into a pair of right-handed neutrinos was the predominant sneutrino 
annihilation channel. In particular, we take the NMSSM parameters as in case C) of 
Table 1 with A Xn = —250 GeV. This correspond to the example studied in Figs. 12 and 
13. In this example the mass of the scalar Higgses is rather large and therefore the 
resulting scattering cross section is reduced. The predictions are plotted in Fig. 19 as a 
function of the sneutrino mass and, as we can see, they are of order 2 x 10~ 9 pb. Once 
more, this result is not excluded by present dark matter searches and can be within 
the reach of some of the future experiments. 

4.1 Overview and comparison with neutralino dark matter 

We have so far shown that the right-handed sneutrino can be a viable WIMP in this 
construction. The examples throughout this Section illustrate how sneutrinos within 
a mass range of xs 5 — 200 GeV and satisfying the WMAP constraint can be 



37 




Figure 19: The same as in Fig. 15 but for example C) of Table 1. 

obtained for which the predicted scattering cross section is not excluded by current 
experimental searches (unlike in the case of the left-handed sneutrino in the MSSM). 
Interestingly, the predicted is generally within the reach of future experiments, 
except along the Higgs resonances in which case the predicted scattering cross section 
can even be below the sensitivity of Ton-scale detectors. 

We can combine the results of this Section in order to determine the region in the 
((7^ , to^J plane that the right-handed sneutrino can cover. This is helpful in order 
to understand in which regions could it account for a hypothetical WIMP detection 
in the future. We do this in Fig. 20, which is constructed gathering all the predictions 
from this Section, and which serves as a summary of our results. 

This also allows us to compare with the predictions for other existing WIMPs. In 
particular, it should be pointed out that the area shown in Fig. 20 is within the area 
predicted for neutralinos in the NMSSM [30, 54] (as well as for neutralinos in the MSSM 
except for the region for very light WIMPs). One may therefore wonder whether these 
two candidates can be discriminated in the case of a hypothetical detection of dark 
matter. It should be pointed out in this respect that since the sneutrino has no spin- 
dependent couplings, the combination of data from two experiments whose targets are 
sensitive to spin-dependent and -independent couplings could be used to disentangle 



38 



1 



10 10 2 

m N1 (GeV) 



Figure 20: Combination of all the theoretical results for the sneutrino elastic scattering 
cross section as a function of the sneutrino mass. 

between sneutrino and neutralino dark matter in the same way that it was shown for 
Kaluza-Klein dark matter and neutralinos in Ref. [55]. 

5 Implication for collider physics 

Finally, we briefly comment on the characteristic signals that might be expected from 
this model at accelerator experiments. 

As with other WIMP models, right-handed sneutrinos can be pair-produced in a 
collider and subsequently escape undetected. This would manifest as events with a 
deficit of transverse momentum. Notice that one expects the same kind of signals if 
the neutralino is the LSP and a dark matter candidate. Thus, the relevant question 
is whether or not the identity of the missing particle (sneutrino or neutralino) can be 
determined. An obvious difference between both particles is their spin. In this sense, 
the study of certain kinematic variables can give us some insight not only about the 
mass of the LSP [56], but even about its spin [57]. This would also be similar to the way 
the nature of the NLSP can be determined in scenarios of gravitino dark matter [58]. 
This, however might not be possible in the Large Hadron Collider and one would need 



39 



to wait for the International Linear Collider. Another crucial difference between the 
sneutrino and the neutralino is their lepton number. This was exploited in scenarios 
of sneutrino NLSP with gravitino LSP by studying the associated lepton production 
[70, 71]. However, this is not directly applicable to our case, since we are dealing with 
a pure right-handed sneutrino. 

Nevertheless, there exists yet another attractive possibility, the production of right- 
handed neutrinos. These can be produced by the decay of either a Higgs boson or a 
neutralino (which can be the NLSP and decay into a sneutrino LSP and a right-handed 
neutrino). Once the right-handed neutrino is produced, it would be long-lived and 
decay, through the tiny mixing with the left-handed component, into a W boson and 
a charged lepton. This would give rise to a displaced vertex which could be observed. 

As we have already stressed, sneutrino dark mater is very sensitive to the properties 
of the Higgs sector. When a Higgs boson with significant singlet component decays, 
depending on its mass, a pair of LSP right-handed sneutrinos and a pair of right- 
handed neutrinos, in other words the missing momentums and long-lived particles, are 
simultaneously produced. This is also a possible unique signal. 

There is another interesting consequence, namely the invisible decay of the Higgs 
boson, since the right-handed sneutrino interacts with the CP-even Higgs through the 
coupling C H o 0i> . This is a common feature of dark matter particles which predominantly 
annihilate through s-channel Higgs exchange and can be detected through t-channel 
Higgs exchange [72, 73]. 

In addition, through collider studies, it is possible to discriminate our model from 
other models with thermal right-handed sneutrino dark matter. In particular, there 
exists a family of models in which thermal right-handed neutralino dark matter is 
viable thanks to the inclusion of a new gauge interaction through a U(l)' [19] or 
U(1)b-l [20]. These scenarios present a Z' gauge boson which is not present in our 
construction. The detection or not of this Z' would discriminate among these models. 
Also, another construction in which the Higgs sector is extended was presented in 
Ref. [23]. The essential difference with our scenario is that the singlet Higgs in that 
model [21, 23] is quite heavy and exchanging Higgs particles are SU(2) doublet only. As 
a result, the correct relic abundance can only be realized with the help of the resonant 
annihilation. Hence, a strong correlation between the sneutrino mass and the Higgs 
mass is present. In addition, the right-handed neutrino production is only through the 
neutralino decay because of too heavy singlet Higgs. Thus, by examining the way of 
right-handed neutrino production, we may discriminate between these two models. 
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6 Conclusions 



We have investigated the properties of the right-handed sneutrino and its viability 
as a WIMP dark matter candidate in an extended version of the NMSSM in which 
a right-handed neutrino superfield is included with a coupling to the scalar Higgs in 
order to provide a Majorana right-handed neutrino mass. Both the \i parameter and 
the neutrino Majorana mass are generated by the vacuum expectation value of the 
singlet field after radiative electroweak symmetry breaking and are therefore naturally 
of order of the electroweak scale. This implies small values for the neutrino Yukawas, 
of order of the electron Yukawa coupling and therefore leads to negligible mixing terms 
in the sneutrino mass matrix and to pure left- and right-handed mass eigenstates. 
The sneutrino spectrum is then analysed in terms of the three new parameters, the 
sneutrino-Higgs coupling \ N , the sneutrino soft mass and the new trilinear term 

We compute the thermal relic abundance for the right-handed sneutrino, including 
all possible annihilation channels. The annihilation cross section is extremely depen- 
dent on the Higgs sector of the NMSSM. Interestingly, light CP-even and light CP- 
odd Higgses are viable in this model as in the NMSSM (if they have a large singlet 
component) and this provides interesting properties for the sneutrino. In particular, 
annihilation into light Higgses (scalar or pseudoscalar) makes it possible to reproduce 
the WMAP result for the relic abundance for low values of the sneutrino mass. We 
study the relevance of each individual annihilation channel and show how sneutrinos 
in the mass range 5-200 GeV can reproduce the correct value for the relic abundance 
for natural values of the three new parameters, = 0(0.1), = (9(100) GeV and 
\A Xn \ = 0(250- 500) GeV. 

Next, the sneutrino-proton elastic scattering cross section is computed and com- 
pared with present experimental sensitivities. We find that the sneutrino is not ruled 
out by current searches for dark matter. Interestingly, many regions of the parameter 
space could be within the reach of some of the future dark matter experiments. 

Finally, we commented on the possible implications for collider physics, identifying 
some potentially interesting signals. We discuss how, in some cases, our construction 
can be distinguished experimentally from other models for sneutrino dark matter as 
well as from other dark matter candidates such as the neutralino. 
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A Feynman rules 



In this appendix we indicate the relevant vertices for the calculation of the sneutrino 
annihilation cross section. The couplings relevant to our computation and which receive 
contributions from the right-handed sneutrino sector we have introduced read 
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also using the following usual NMSSM Feynman diagrams. 
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The explicit expressions for the couplings C H o H o H o C a oao h o C h + h - h o C h o a o z , 

i j k a, b k k fee, 

and C H o H+w - can be found in [59], and translated into our formalism by doing the 
transformations v\p — > ^= and v s — > 
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B Annihilation channels 



In order to determine the total annihilation cross section we follow the procedure of 
[38] and start by introducing a Lorentz-invariant function w(s) 



w(s) 



dUPS\A(NN ^all)| ; 



(B.29) 



where \A(NN — > all) | 2 denotes the absolute square of the reduced matrix element for 
the annihilation of two sneutrinos, averaged over initial spins and summed over final 
spins. The function w(s) is related to the annihilation cross section a(s) via [61] 



w(s) = ^s(s-Am^ L 2 )a(s). 



(B.30) 



Since w(s) receives contributions from all the kinematically allowed annihilation 
process NX — > XiX 2 , it can be written as 



wis 



) = — £ 

32 7T ^ 



X\ X 2 



c9 (s - (m Xl +m X2 ) 2 ) /3 x (s,m Xl ,m X2 )w XlX2 (s) 



,(B.31) 



where the summation extends over all possible two-body final states XiX 2 , m Xl and 
m X2 denote their respective masses, and 



c = 



c f \{X m =f{j) 
1 otherwise, 



^B.32) 



where C/ is the color factor of SM fermions (cj = 3 for quarks and Cf = 1 for leptons). 
The kinematic factor (5 X is defined as 



p x (s,m Xl ,m X2 ) 



(m Xl +m X2 f 
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^B.33) 



In the CM frame, which we choose for convenience, the function w XlX . 2 (s) can be 
expressed as 



w XlX2 {s) = 



+i 



d cos 6 CM \A(N X -> X!X 2 )| 2 , 



^B.34) 



where 6 C m denotes the scattering angle in the CM frame. In other words, we write 
\A(XX — > X!X 2 )| 2 as a function of s and cos9 C m, which greatly simplifies the com- 
putation. 
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We present here the analytic expressions for the w functions for the various anni- 
hilation channels. In doing so, we have assumed the lightest sneutrino to be a pure 
right-handed sneutrino, hence Nf~ = and NY^ = 1 with % — 1. 

In the following expressions we define 



, ('/>//< I\) 2 



mjjo) + , (B.35) 



s - m H0 



= (s — m 2 H o)(s — m 2 H o) + m H o m H oTiTj 



(m H o Ti(s - m 2 ) - m H o TJs - m 2 )) 2 
+ ' 2 w- ' , ,1 - ^ , (B.36) 



^ p )(s - «^p) + m H° FiTj 



for i, j = 1, 2, 3, where Tj is the decay width of the Higgs iff. 



The annihilation proceeds through an s-channel mediated by a scalar Higgs, iff. 

ww+w- = Qg 2 M w 2 — ±7. ~ (B.37) 

i,j=i v 



• ^ zz 

The annihilation proceeds through an s-channel mediated by a scalar Higgs, iff. 
The t-channels where a sneutrino is exchanged do not contribute since the Z — v — v 
coupling vanishes for pure right-handed sneutrinos. 

Sg 2 M z 2 A C H?ui> C H o^C Hi C Hj 

w Z z = ^— > j. B.38 

cos9 w fp Aij 



There are four contributions to this annihilation channel. One comes from an 
s-channel mediated by a scalar Higgs, iff, the t, w-channels where a sneutrino is ex- 
changed, and the point interaction. 



48 



2{Lh°h°vv) - 4 2^ a + 2 » + 



fe=l 



* J i i 

A^A^E^s 

i i 1 

( Tit 2 . jDU 2\ ( TDt TDU \2 

(l-B^il-B^iB^+B^: 



k,l=l 

^i(-B^-p) + F\{B U H0 



ki 



r>t r>u 



+ 



C h° H° H°i>i>C H q pC fjO^A^o + A u „ ) 



fc=l 



i i % 

C H?vi>C H°.vi>C H"vv(.A t „ + -A" ) 

) * £ z i i 

' .v lv . r;, ,/•;„= v* 



^l(^ )+^ 1 (^o) 
i g 



^i(-Bhp) + FiiB^o] 

4 Z_ 



(B.39) 



which should be multiplied by a factor 1/2 in the case of identical particles in the final 
state. 

We have defined 



^(x) = log 



(B.40) 



and 



where 



<o = -7=^-1 (B.41) 



m 2 H o - s 



A U H0 ee -7^— + 1 (B-42) 
\fsh H o 



Vi- m Ni 2 \/ E H? 2 - m ] 

yt — V 



£^0 = — 2 ' ~ ( B - 43 ) 

m^ - y/sEjjO 



H? 



B > = 2 " 7=73 L (B.44) 
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•s - (m 2 H o - m 2 H ^ 



' ( m 2 H0 



The annihilation proceeds through an s-channel mediated by a scalar Higgs, H® 
and a point interaction. 



W A o A o = 2 { C A°A°i>i>) - 4 2^ A + 

k=l k 

2 ^ C HpuCH^uuC A o A O bH oC A o A o bH o ^ ^ 

k,l=l 

which is multiplied by a factor 1/2 in the case of identical particles in the final state. 



•->// 

We only include here the s-channel mediated by a scalar Higgs, H®. Once more, 
the s-channel annihilation mediated by the Z boson vanishes for a pure right-handed 
sneutrino. For annihilation into 11 (uu) one should, in principle, also consider the t- 
channel mediated by charginos (neutralinos). However, those diagrams are suppressed 
by the very small Yukawa coupling, Yn, and will be here neglected. 

— : \ - 3 (2a - 8m}) (B.48) 



• -> NN 

There are three contributions to this channel, namely the s-channel, mediated by 
a scalar Higgs, H®, together with the t- and u- channels where (the five) neutralinos 
are exchanged. 
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First, we address the case where the lightest sneutrino is the real part of right- 
handed sneutrino. 

s i st,su . t.u /t-) Ar.\ 

w NN = w NN + w N ' N + w^ N , (B.49) 

with 



w s 



ST C HpuCH«vuC NN HlC NNH » AM 2s 

a \ - 2^ a y s ~~ M ' V ) 

k,l=l 



(B.50) 



w 



nn - 2 / 2 / 



fc=l 1=1 



P(?7lc.) , 



A; - A J 



(B.51) 



ID 



NN 



32C ON j,. 



~ + 



i=l 

5 /nr2 

E°^ONxi iWjt 
S 2 



s 2 t A(m^) B{m jii ),(A i -A\ 



8 (A 2 - A 2 ) 2A { A 
s 2 C(m^,m^) 



log 



Ai + A 



+ 



log [^4 1 



A (A 2 -A 2 ) \Ai - A 



log ' 



(B.52) 



The following auxiliary functions have been used 
A(m^) = -m^) (m N + m^f + | -m^ 2 + 2m 2 N A i + 2m N m^A^j - 

fjU- > (B.53) 
B(m i4 ) = (| -m^) (mjv + m^) 2 - | (| - m^ 2 + 2m 2 v A i + 2mjvm^ + 

3^4 2 , (B.54) 
C(m^,m^.) = (| -m%) (m N + m^) (m N + m jij )A j + 
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We have also defined 



A = -^s-A mfl 2 ^js-Am 2 N , 



(B.55) 
(B.56) 

(B.57) 
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In the case that the lightest sneutrino is the imaginary part of right-handed sneu- 
trino, one can obtain the relevant w^n by replacing 

st.su st.su /t-j rrv\ 

W NN ~ W NN , ( B - 59 ) 

and 



m N — > — m N (B.60) 



m w 



t,u 
NN- 



• -> H+H- 

As in the case of annihilation into pseudoscalar Higgses, the annihilation receives 
contributions from an s-channel mediated by a scalar Higgs, and a point interaction. 



2 A Ch+h-ppChOmCh+H-HQ 
WH+H- = ^H+H-vu) - 4 > r h 

k=l K 



This channel receives contribution from s-channel scalar Higgs exchange. Notice 
that although there would also be an s-channel with Z boson mediation, the smallness 
of the Yukawa Y N renders that contribution negligible. 



WZAO 



3 fi r~t r~i /~i 



k,l=l 



s 2 + m mi 4 — 2m ao 2 s\ 
-2s + M z 2 + 2m A o 2 + ^2 — . (B.62) 



M 



z 
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• -> W~H + 

As in the case above, this channel receives contribution from s-channel scalar Higgs 
exchange. 



EC Hlvi>C H°Oi>C H°H+ W~ Ch°H+W~ 




s 2 + m H + 4: — 2m H +' 2 s 
M w 



2s + My/ + 2m H+ 2 + : g • ( B - 63 ) 
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